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An exact correspondence is established between a A^-body classical interacting system and a — 1- 
body quantum system with respect to the partition function. The resulting quantum-potential is a 
N — 1-body one. Inversely the Kelbg potential is reproduced which describes quantum systems at a 
quasi-classical level. The found correspondence between classical and quantum systems allows also 
to approximate dense classical many body systems by lower order quantum perturbation theory 
replacing Planck's constant properly by temperature and density dependent expressions. As an 
example the dynamical behaviour of an one - component plasma is well reproduced concerning the 
formation of correlation energy after a disturbance utilising solely the analytical quantum - Born 
result for dense degenerated Fermi systems. As a practical guide the quantum - Bruckner parameter 
rs has been replaced by the classical plasma parameter F as « O.SF^''^. 



Several hints in recent literature conjecture that there 
seem to exist a correspondence between quantum systems 
and higher dimensional classical systems. The authors of 
[01 argue that a higher dimensional classical non-Abelian 
gauge theory leads to a lower dimensional quantum field 
theory in the sense of chaotic quantisation. The corre- 
spondence has been achieved by equating the tempera- 
ture characterising chaotization of the higher dimensional 
system with h of the lower dimensional system by 

h = aT. (1) 

Recalling imaginary time evolution as a method to cal- 
culate correlated systems in equilibrium such correspon- 
dence seems suggestible. We will find a similar relation 
as a best fit of quantum - Born calculations to dense in- 
teracting classical systems. 

In condensed matter physics it is a commonly used 
trick to map a two - dimensional classical spin system 
onto a one - dimensional quantum system This sug- 
gests that there might exist a general relation between 
classical and higher dimensional quantum systems. We 
will show that a classical many body system can be 
equally described by a quantum system with one particle 
less in the system but with the price of complicated non- 
local potential. This can be considered analogously to the 
Bohm interpretation of quantum mechanics where the 
Schroedinger equation is rewritten in a Hamilton- Jacobi 
equation but with a nonlocal quantum potential. 

Another hint towards a correspondence between classi- 
cal and quantum systems was found recently in Q where 
it was achieved to define a Lyapunov exponent in quan- 
tum mechanics by employing the marginal distribution 
which is a representation of Wigner function in a higher 
dimensional space. Since the Lyapunov exponent is es- 
sentially a concept borrowed from classical physics this 
finding points also in the direction that there exists a 
correspondence between quantum systems and higher di- 
mensional classical systems. 



On the opposite side there are systematic derivations 
of constructing effective classical potentials such that the 
many body quantum system is described by the classical 
system. An example is the Kelbg potential for Coulomb 
systems [0-|| 

F--^W.^(l-e-^V.VV^^erfc(9) (2) 

with P = h? /2iiT and = l/mi-|-l/m2 describing the 
two-particle quantum Slater sum correctly by a classical 
system. Improvements and systematic applications can 
be found in 

Here in this paper it should be shown that a classical 
A^-particle system can be mapped exactly on a quantum 
N — 1-particle system with respect to the partition func- 
tion. Though the resulting effective iV — 1 body quantum 
potential is highly complex it can lead to practical ap- 
plications for approximating strongly correlated classical 
systems. In the thermodynamical limit it means that the 
dense classical system can be described alternatively by 
a quantum system with properly chosen potential. 

This finding suggests that the quantum calculation in 
lowest order perturbation might be suitable to derive 
good approximations for the dense classical system. This 
is also motivated by an intuitive picture. Assume we 
have a dense interacting classical plasma system. Then 
the correlations will restrict the possible phase space for 
travelling of one particle considerably like in dense Fermi 
systems at low temperatures where the Pauli exclusion 
principle restrict the phase space for scattering. There- 
fore we might be able to describe a dense interacting 
classical system by a perturbative quantum calculation 
when properly replacing fi by density and temperature 
expressions. Indeed we will demonstrate in a one - com- 
ponent plasma system that even the time evolution and 
dynamics of a very strongly correlated classical system 
can be properly approximated by quantum - Born cal- 
culations replacing the quantum parameters by proper 
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classical ones. 

Let us now start to derive the equivalence between clas- 
sical and quantum systems by rewriting the classical N- 
particle partition function. The configuration integral 
reads 



N 

Qn{P) = / dxi...dxN ]^(1 + fij) 



(3) 



where we used 

Meyer's graphs fij = cxp {—f3uij{xi — Xj)) ~ 1 with the 
interaction potential (xi — Xj ) of the classical particles 
and the inverse temperature /3. It is now of advantage to 
consider the modified configuration integral 

Qn{P)^Qn{2P) 



(4) 



X (1 + /l2)(l + /l3)(l + /l4)...(l + flN) 

X (1 + /21')(1 + /23)(1 + /24)...(1 + f2N) 



X (1 + fNv){l + /JV2')(1 + /W3')-(1 + fNN-v) 

such that a quadratic schema in fij' appears. Now we 
assume a complete set of — 1 particle wave functions 
^njv_i such that 

5{xi - x[)...S{xN-i - ^^'at-i) 
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with some "quantum numbers" {i} characterising the 
state. Further we propose the following eigenvalue prob- 
lem defining the wave function 



N 

/ dxi ]^(1 + /ij)^'ii..ijv_i(xi...XAr_i) 



(6) 



with the system volume V. This allows to calculate the 
configurational integral (^) exactly by successively inte- 
grating Xi.-.Xn 
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= Yl j dxi...dxNdx[...dx'N-i-i'*,..^„_,ix[...x'j^_i) 

■ii..ijv-l 

'X-'^n..iN-iixi---XN-l) 

X (1 + /l2)(l + /l3)(l + /l4)...(l + flN) 
X (1 + /21')(1 + /23)(1 + /24)...(1 + f2N) 
X 

X (1 + /Ari,)(l + /iV2')(l + /JV3')-(1 + fNN-v) 



xy^e-^-<->*.,...„_,(x;...:r^_i) 
yN J2 e-^^{'>. 



(7) 



This establishes already the complete proof that we can 
map a classical TV-body system on a A'^— 1-body quantum 
system since (^) is the eigenvalue problem of a — 1- 
body Schrocdinger equation. To see this we can consider 
a wavefunction ^ built from the Fouriertransform of 'J 
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_^ipi...pN~i) (8) 



which obeys the N — 1-particle Schrocdinger equation 




with Efq^i cx Ve ^f'> and we rewrote the left hand side 
of (j^) as quantum potential 



< xiii...XN-iiN-i\U\x\i\...Xf^_ii'f^_i > 

^I3[ui2(x[-Xi) + ...+Uin{x[-Xn-i)] 

xn6{xi - x'2)...6{xn_2 - (10) 



The resulting equivalent quantum potential (|1C|) is a 
N — 1-body nonlocal potential with respect to the co- 
ordinates but depends on N strength function parame- 
ter (e.g. charges). Therefore we have casted a classical 
A^-body problem into a nonlocal quantum A^ — 1 body 
problem. One could easily give also a symmetrised or 
anti-symmetrised form of the potential using symmetries 
of the wave function and permuting coordinates of ( |lO| ) 
respectively. We do not need it here since we will restrict 
to applications neglecting exchange correlations further 
on. 

While the above correspondence holds for any particle 
number and might be useful to find solvable models for 
classical three - body problems, we will consider in the 
following many - body systems. First let us invert the 
problem and search for an effective classical potential ap- 
proximating quantum systems. This should us lead to the 
known Kelbg-potential (||). For this purpose we assume 
a quantum system described in lowest approximation by 
a Slater determinant or a complete factorisation of the 
many - body wave function into single wave function 
^h...iN{xi...XN) = (t)ii...(t>i^. We neglect for simplicity 
exchange correlations in the following. The correspond- 
ing eigenvalue equation for (j) itself one can obtain from 
(I) or d) by multiplying with ^f*^ j^_^(x2...XAr_i) and 
integrate over X2.--Xn-i- To see the generic structure 
more clearly we better calculate the correlation energy 
by muhiplying (|) or (|) by ^f*^ ^^_^,^(x2...a;jv) and in- 
tegrating over X2---XN- This provides also the eigenvalue 
e{i} and leads easily to approximations for the partition 
function (|^). To demonstrate this we choose the lowest 
order approximation taking identical plane waves for <j). 
Than the pressure can be obtained from the partition 
function via (0) 
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P-T ^ hnO N{N~2) 



(11) 



where V is the volume of the system. We recognise 
the standard second virial coefficient for small potentials 
while for higher order potential the factor 1 /2 appears in 
the exponent instead as a pre-factor indicating a different 
partial summation of diagrams due to the schema behind 
(0) and d). 

To go beyond the plane wave approximation we multi- 
ply (||) by ^.*2-.jjv-iji (^2---a;Ar) and the kinetic part of the 
statistical operator before integrating over X2---Xn- This 
means we create an integral over the A'^ — 1 particle den- 
sity operator and the potential ( [lO| ) which together repre- 
sents the correlation energy. This expression is a succes- 
sive convolution between the cluster graphs fij and the 
relative two - particle correlation function p^^^a (^i ~ ^2)- 
The resulting mean correlation energy density reads 

U MT^ f dyi..dyN-i r \ r \ 1 \ 



x(l-/l2(yi))(l-/l3(yi+2/2))-.-(l-/lJv(yi + --- + 2/Ar-l)) 

f dyi..dyN- 



yN- 



xui2{yi)ui3{yi+y2)-.uiNiyi + ---+yN-i) + ■■ 



(12) 



in dimensionless units where all other cluster expansion 
terms lead either to lower mean field or disconnected 
terms. While these terms can be calculated as well we re- 
strict to the highest order convolutions in the correlation 
energy (|l^) which have now the structure of mean corre- 
lation energy U /V — ^ J ^Pixi^ (2^)^12^ with a classical 
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effective potential V^f^ 



3 



Vf{r) (X 



34 



V 

dxidx2 



Pl2{xi)ui2{xi)u23,{xi + r) (13) 
Pl2{xi)ui2{xi)uiz{xi +X2) 



1/2 

xp23(a;2)u34(a;i + X2 + r) 



(14) 



where the two-particle, three-particle etc. approximation 
can be given. In equilibrium the nondegenerate correla- 
tion function reads [P — h? / ^iT = /2ii\ 



Pilt2 ixi-X2)'' 



dp 



(15) 



Using the Coulomb potential m cx l/r we obtain from 
the two-particle approximation ( |l3|) just the Kelbg po- 
tential (H). The three - particle approximation (|^ can 
be calculated as well and reads [x — r/l] 



(16) 



The comparison of the third order potential with the 
Kelbg potential can be seen in figure |^. The third order 
potential is somewhat less bound than the Kelbg poten- 
tial. 




FIG. 1. The comparison of the Kelbg potential (g) and the 
third order potential ( p^ versus r/l. 

With the schema (|l^ one can easily integrate higher 
order approximations as successive convolutions, but 
with respect to the small differences between (|^) and ( p^ ) 
in figure ^ one does not expect much change. Also in 
principle the degenerate case could be calculated using 
Fermi-Dirac distributions in ([l5|). But one should then 
consider also the neglected exchange correlations during 
factorisation of 5" as well. Let us summarise that the 
known effective classical potential describing a quantum 
system in binary approximation has been recovered by 
identifying the effective two-particle interaction within 
the correlation energy. 

We want now to proceed to a phenomenological level 
in that the above correspondence between quantum and 
classical systems motivates to find good approximations 
for the dynamics of classical many-body systems by em- 
ploying quantum-Born approximations. This can be un- 
derstood by the fact that the Kelbg potential deviates 
appreciably from the Coulomb one only if the interparti- 
cle distance d are smaller than the thermal wave length 
A. In other words for dense classical systems under such 
conditions we can think of it as a dilute quantum system 
replacing A ~ d. To check this conjecture let us consider 
an one-component plasma system which is characterised 

by two values. The classical coupling is described by 

2 

the plasma parameter F = as a ratio of the length 
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where Coulomb energy becomes larger than kinetic en- 

2 

ergy, to the interparticle distance or Wigner size ra- 
dius d = Ideal plasmas are found for F << 1 
while around F = 1 non- ideal effects become important. 
A second parameter which controls the quantum features 
is the Bruckner parameter as the ratio of the Wigner size 
radius to the Bohr radius as — /me^. Quantum ef- 
fects will play a role if < 1. We will consider the 
situation that the interaction of such system is switched 
on at initial time. Then the correlations are formed by 
the system which is seen in an increase of temperature 
accompanied by the build up of negative correlation en- 
ergy. This theoretical experiment has been investigated 
numerically by ||l2| for classical plasmas with different 
plasma parameter F. 

In 1^,0 we have calculated the formation of such 
correlations by using quantum kinetic equations in Born 
approximation. The time dependence of kinetic energy 
was found at short times to be 



dkdpdq 



1 



ab ^ ' 

y<f'anil-fa){l-fb) 



(17) 
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where / are the initial distributions and 
J) (p+q) ^j^g statical screened Coulomb in- 

2mfa 2ma Inib 

teraction is Voiq) — 'ine^h^ / {q^ ~\-h^K'^) with the inverse 
screening length expressed by density n and temperature 
T as = Aire'^n/T or — Qire^n/ef for the high or 
low temperature limit. For both cases dynamical as well 
as statical screening it was possible to integrate analyt- 
ically the time dependent correlation energy ([l7|). This 
has allowed to describe the time dependence of simula- 
tions in the weak coupling limit F < 1 appropriately [p^ . 
For stronger coupling F > 1 the Born approximation fails 
since the exact correlation energy of simulation is lower 
than the first order (Born) result kc^ /2T = ^JJj2T^/'^. 
Moreover there appear typical oscillations as seen in fig- 
ure ||. 

Now we will employ the ideas developed above and will 
use the quantum Born approximations in the strongly 
degenerated case to describe the classical strongly corre- 
lated system. For strongly degenerated plasmas the time 
dependence of correlation energy was possible to inte- 
grate as well with the result expressed here in terms 
of plasma parameter F and quantum Bruckner parameter 



nT 



sm yr 



(367r4)i/6 T \ yr 

X f larctan(l' 
\oi bi' 



bf + bf 



(18) 



with bi = hK/2pf = Vr/(487r2)i/6^ yj- = Atft/n = 
(2)''/^7r^/^3^/^r/Y^ where the time is scaled in plasma 



periods t = 27rt/ajp. Now we try to fit this quantum 
result to the simulation using the Bruckner parameter 
as free parameter. For the available simulations between 
1 < F < 10 we obtain a best fit 




0.5. 



(19) 



The quality of this fit is illustrated in figure which is 
throughout the range 1 < F < 10. This is quite astonish- 
ing since not only the correct classical correlation energy 
|l5| is described but also the correct time dependence i.e. 
dynamics. 




t [27l/(Op] 

FIG. 2. The time evolution of a classical one-component 
plasma after sudden switching of interaction jl^] compared 
to the quantum Born result when the Bruckner parameter is 
replaced according to (^^. The long time equilibrium value is 
remarkably well reproduced by the quantum - Born result (p^. 

Let us try to understand what this phenomenological 
finding means. Using the thermal De Broglie wave length 
= /AmT we can rewrite (fol) as 



A^ 



d 

K 



(3F)i/2- 



(20) 



In the considered range of F — 1...10 we have (3F)^/^ = 
1...2 and the thermal wave length A is found to be nearly 
equal the interparticle distance d as a best fit of quantum 
Born calculation to dense classical systems. This is ex- 
actly the distance where the Kelbg potential ( p^ or (||) 
starts to deviate from the Coulomb potential. In other 
words we confirm the conjecture that the dense classical 
system can be described by dilute quantum systems if in 
the latter systems the thermal wave length is replaced by 
the interparticle distance. This condition (^9|) can also 
be rewritten into the result (^ of literature using the 
degenerated screening length. 

We summarise that in equilibrium we have shown that 
there exist an exact relation between a A^-body classical 
system and a N — 1-body quantum system. This has 
allowed to recover the quantum Kelbg potential easily. 
As practical consequence we suggest to describe the dy- 
namics of dense interacting classical many body systems 
by the simpler perturbative quantum calculation in de- 
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generate limit replacing properly Ti by typical classical 
parameters of the system. 

I would like to thank S. G. Chung for numerous dis- 
cussions and valuable hints. 
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